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ABSTRACT 


We  give  a treatment  of 
molecule-ion  molecule  reactions 
based  on  stochastic  mechanics. 
Stochastic  mechanics  is  a semi- 
classical  theory  in  which  one 
assumes  that  particles  move  as  a 
stochastic  diffusion  process.  It 
is  shown  that  the  reaction  prob- 
ability and  reaction  rate  can  be 
determined  from  the  solution  of 
certain  partial  differential  equa- 
tions (pde) . Asymptotic  solutions 
(for  h-»-0)  of  these  pde  are  con- 
structed in  terms  of  incomplete 
special  functions.  The  results 
derived  using  stochastic  mechanics 
are  compared  with  results  derived 
using  other  semiclassical  approxi- 
mations . 
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SECTION  1. 


INTRODUCTION 


In  this  paper  we  consider  the  reaction  of  molecule  ions  and 
molecules  from  the  viewpoint  of  stochastic  mechanics  (1) . The 
reactions  of  interest  are  of  the  form 

+ YB  - XB'*’  + Y,  (1.1) 

e . g . 

Ar'*'  + -*  ArH"^  + H.  (1.2) 

These  reactions  may  lead  to  the  phenomenon  or  orbiting  (see 
below)  and  have  been  studied  by  Hirschfelder  et.  al.  (2) , Gioumousis 
and  Stevenson  (3)  and  Ford  and  Wheeler  (4).  These  wor)cers 
analyzed  the  reaction  (1.1)  by  using  classical  or  semiclassical 
]cinetic  theory  (2,3)  or  by  taking  the  semiclassical  limit  of  the 
quantum  mechanical  scattering  theory  (4) . The  object  of  major  im- 
portance in  these  theories  is  the  reaction  cross-section,  o.  once 
the  cross-section  is  known,  reaction  rates  and  probabilities  can  be 
calculated.  In  (4),  rainbow  and  glory  scattering  are  given  rela- 
tively simple  and  general  semiclassical  treatments.  The  authors 
conclude  that  "there  is  no  semiclassical  approximation  to  the 
orbiting  effect  of  simplicity  or  generality  comparable  to  the 
analysis  for  rainbow  scattering  and  glory  scattering"  (4,  page  274). 
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In  this  paper,  we  give  a general,  and  intuitive,  semi- 
classical  approximation  to  orbiting.  Our  approach  is  very  dif- 
ferent from  the  approach  in  (2-4)  , in  that  stochastic  mechanics 
(SM)  is  used.  One  advantage  of  the  SM  approach  is  that  reaction 
rates  and  probabilities  can  be  calculated  directly,  once  the 
stochastic  problem  is  formulated. 

In  section  2,  the  classical  mechanics  problem  is  discussed. 

A discussion  of  the  phenomenon  of  "orbiting"  is  given.  In  section 
3,  the  stochastic  mechanics  problem  is  formulated.  Stochastic 
mechanics  assumes  that  a particle's  motion  is  a diffusion  process. 
The  drift  term  of  the  diffusion  is  given  by  classical  mechanics 
and  the  intensity  of  the  diffusion  is  fi/2m,  where  m is  the  mass  of 
the  particle.  The  reaction  probability  and  reaction  rate  are  then 
determined  by  solving  certain  elliptic  partial  differential  equa- 
tions. The  solution  of  these  equations  is  discussed  in  section  4. 


SECTION  2 


CLASSICAL  MECHANICS  OF  MOLECULE  ION-MOLECULE  REACTIONS 

We  assume  that  the  YB  molecule  is  chosen  as  the  origin. 

The  molecules  approach  the  YB  molecules  with  initial  velo- 

city g and  asymptotic  center-to-center  distance  b (see  Fig.  1) . 
The  center-to-center  distance  b is  called  the  impact  parameter. 

We  describe  the  collision  process  by  the  use  of  r(t),e(t)  co- 
ordinates, shown  in  Fig.  1.  According  to  classical  scattering 
theory  (2,3,5)  r(t)  and  0 (t)  evolve  according  to 


dE  g 


A - "9^ 

r 


r 


Zae 

2r^E 


(2.1) 


(2.2) 


where  Z is  the  valence  of  X , a is  the  polarizability  of  YB 

and  E is  the  initial  energy  of  X^.  The  molecule  ion  will  "orbit" 
• • 
when  r = 0.  The  condition  for  r = 0 is 


2 2 
, b^  ^ Zae 

1 7 — 4“ 

r 2r  E 


= 0. 


(2.3) 


This  equation  will  have  a real  solution  only  if 


b > b = J 2Zae 
— c 


2 \ 1/4 


(2.4) 
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We  call  the  critical  value  of  the  impact  parameter.  If 

b < b^,  then  all  solutions  of  (2.1,  2.2)  will  spiral  into  the 
origin.  Once  the  molecule  ion  is  within  a criticax  reaction  dis- 
tance, r^,  of  YB,  we  assume  that  a reaction  occurs  with 
probability  1.  If  b ^ b^ , then  more  interesting  behavior  may 
occur.  For  b > b^,  the  molecule  ion  will  approach  YB  and  then 
fly  away  on  a hyperbolic  orbit.  When  b = b^ , however,  the  mole- 
cule ion  will  be  trapped  in  an  orbit  around  YB.  The  radius  of 
this  orbit  is  calculated  from  (2.3)  and  is 


We  view  (2.1,  2.2)  as  a classical  dynamical  system.  The  orbit  at 
r = r^  is  an  unstable  limit  cycle  (see  e.g.  (6)  for  a discussion 
of  these  terms) . 

In  any  experiment,  the  molecule  ion  will  be  produced 
according  to  some  initial  distribution  of  impact  parameters, 
pQ(b).  Using  the  classical  analysis,  we  see  that  the  reaction 
probability  will  be 


The  time  that  it  takes  to  enter  the  reaction  sphere  r^  from  a 
distance  r^  is 


T(b) 


1 


g dr 

.2  „ 2-11/2 
b , Zae 

r 2Er 


(2.7) 


Hence,  the  mean  time  to  enter  the  reaction  sphere,  averaged 
over  those  molecules  that  enter  is 


<T(b)  > 


T(b)  Pq  (b)db 


(2.8) 


The  classical  reaction  rate  is  defined  to  be 


k 


1 

<T(b)>  • 


(2.9) 


If  the  system  involved  more  complex  potentials,  more  compli- 
cated dynamical  systems  would  be  derived.  For  example,  for  the 
Lennard-Jones  potential 


V(r)  = 4n 


equation  (2.3)  is  replaced  by 


(2.10) 


(2.11) 
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In  this  case,  stable  and  unstable  limit  cycles  will  appear  in  the 
classical  dynamics  (Fig.  2). 
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SECTION  3. 


STOCHASTIC  MECHANICS 


In  order  to  convert  the  classical  theory  of  section  2 to  a 
semiclassical  theory,  we  assume  that  r(t)  and  ^I’Ct)  are 
stochastic  diffusion  processes  (1).  Specifically  assume  that 
f(t)  and  ^(t)  satisfy  the  Ito  equations 


In  (3.1,  2)  and  are  independent  Brownian  motion  processes. 

2 

Namely  E(dW)  = 0,  E(dW  ) = dt  and  the  increments  dW  are 
normally  distributed. 

The  picture  provided  by  stochastic  mechanics  is  this:  par- 
ticles move  according  to  classical  mechanics  and  are  continually 
perturbed  by  quantum  noise. 

Let  u(r,0),  T(r,0)  be  defined  by 


u{r,e)  = 


T(r,0)  = E 


Pr  I enters  r^  rather  than  approaching  r = “> 

|r(0)  = r,?(o)  = ej  (3.3) 

< t:  r(t)  £ r , r(s)  > r for  s < t|r(0)  = r 
^(0)  = 0 J ^ (3.4) 
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The  function  u(r,0)  is  the  reaction  probability  as  a function  of 
initial  position  (and  thus  impact  parameter).  T(r,9)  is  the 
average  time  that  it  takes  a particle  to  enter  the  reaction 

sphere  of  YB,  given  the  initial  position.  Hence,  after  an  appro- 
priate averaging  over  initial  position,  T(r,6)  is  a measure  of 
the  rate  of  the  reaction. 

In  a standard  way  (1,7)  it  can  be  shown  that  u and  T 
satisfy  the  following  deterministic  partial  differential  equations: 


jh 

2m 


(3.5) 


jh 

2m 


(T 


rr 


T ) 


-1. 


(3.6) 


In  (3.5,6)  subscripts  indicate  partial  derivatives. 

The  boundary  conditions  for  (3.5)  are 

u(r  ,9)  = 1 u(<»,9)  = 0 u(r,9  + 2tt)  = u(r,9).  (3.7) 

The  boundary  conditions  for  (3.6)  are 
T(r^,e)  = 0 lim  T^(r,9)  = 0 


T(r,9  + 2Tr) 


T(r,9)  . 


(3.8) 


Once  u(r,0),  T(r,e)  are  known,  the  overall  reaction  prob- 
ability and  rate  are  obtained  by  averaging  over  the  initial 
distribution 

V = Ju (r, 0) pQ (r, 0)drde  (3.9) 

k = ^|T(r,0)  Pq  (r,0)drd0j"^  . (3.10) 

Equations  (3.5,  6)  are  singular  perturbation  problems,  since 
fi/2m  is  small. 
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The  error  intergral  satisfies 


d^E 


dE 


dz 


2 dz 


(4.5) 


In  order  to  evaluate  ip,g  and  f,  derivatives  of  (4.1)  are 
evaluated  and  substituted  into  (3.5).  When  evaluating  derivatives, 

(4.5)  is  used  to  replace  E' 

After  substitution  into  (3.5),  terms  are  collected  according  to 


i'li/yj  fi/2m)  by  -^i|^\/(^i/2m)^  E 


powers  of  ('fi/2m)  . The  coefficients  of  the  various  powers  of 
(-fi/2m)  are  set  equal  to  zero,  yielding  equations  for  ip , g and 
f^.  Here  we  only  discuss  i|',g^  and  f^,  which  gives  the  leading 
term  in  (4.1).  Other  terms  are  treated  in  an  analogous  fashion  (8) 
We  set 


b^  = zM 


(4.6) 


Zae 

i 

2Er 


(4.7) 


• 0 • ^ 

so  that  the  classical  equations  are  9 = b , r = b . The  leading 


term  in  the  asymptotic  solution  of  (3.5)  is 
vanishes  if  ij;(r,6)  satifies 


0 ^(^^/2m)"^/^^. 


It 


+ b%Q  . ^ 


(4.8) 
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The  0(1)  term  vanishes  if  g^(r,9)  satisfies 

b^g°  + h\l  = 0.  (4.9) 

The  0 ^(h/2m)  term  vanishes  if  f^(r,6)  satifies: 

» b"t“  * + ♦jg)  + 2j^gj4._.  + (4.10) 

- - £®((*»,)r  * IWele)  = “ 

We  first  consider  (4.8).  The  region  of  major  interest  is  a 
vicinity  of  the  unstable  limit  cycle,  so  that  we  will  obtain 
iplr,0)  in  a vicinity  of  the  limit  cycle.  Since  r^  corresponds  to 
an  orbit,  b (r^)  = 0.  Intuition  suggests  that  we  set 

= 0.  Then  (4.8)  indicates  that  ip  = 0 on  r^.  We  can 
determine  i;j(r,Q)  in  a vicinity  of  r^  by  Taylor  expansion, 
once  i|)^(r^,0)  is  known.  If  (4.8)  is  differentiated  with  respect 
to  r and  evaluated  on  r = r^ , we  obtain 

b'^ij;„  + bf  i|j  - = 0.  (4.11) 

rer  'r  r r 

Now  b®  9/90  + b^  9/9r  = d/dt,  where  t measures  "time"  along 

the  limit  cycle.  Since  b^(r^)  = 0 and  0 (t)  = -(bg/r^)t  + Qg  on 

r = r , we  obtain 
c 

^i|)+bfi(i-iji^  = 0.  (4.12) 

dt  r 'r^r  r 
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i 


In  tern.s  of  6, 


we  obtain 


2 

r 

c 


+ bf. 


r 


(4.13) 


The  periodic  solution  of  (4.13)  gives  'i>^(r^,0).  Once  this 
derivative  is  known,  •p(r,0)  can  be  determined  in  a vicinity  o‘ 
the  limit  cycle.  A more  accurate  procedure  for  the  determinati.  r 
of  p(r,8)  would  involve  the  solution  of  (4.8)  by  the  method  of 
characteristics  (see  (8)).  Since  b is  independent  of  6,  the 
periodic  solution  of  (4.13)  is  simply 


P = (bf  ) 

r r 


1/2 


(4.14) 


Equation  (4.9)  indicates  that  is  constant  on  classical  trajec- 

0 


tones.  We  set 


identically  constant,  and  determine  its  value  b 


by  requiring  g^E  (ijj  (r  , 0 ) / y(^/2m)  asymptotically  satisfy  the  boundary 
conditions.  From  equation  (4.10),  we  can  derive  an  equation  for 
on  the  limit  cycles  requiring  that  f*^(6  + 2^)  = f*^(6).  Once 
f^  is  known,  the  first  term  of  the  asymptotic  solution  of  (3.5) 
is  known.  From  (4.1) 


u(r,e) 


((b^. 


^ 1/2  , 

) (r 


^fi/2mj  . 


(4.15) 
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of  fall  off 
parameters : 

(4.16) 

The  result  (4.16)  is  a SM  analogue  to  equation  (4.3)  of  Ford  and 
Wheeler  (4) . 

The  solution  of  equation  (3.3),  for  T(r,0)  can  be  given  in 
a form  analogous  to  (4.1): 

T(r,e)  = g (r  , 9 ,'fi/2m)  F (r  ,Q)  / ■^'h/ 

+ (4.17) 

+ k (r , 0 ,ii/2m)  . 

In  (4.17),  g,f  and  4)  are  defined  as  above,  k (r , 0 ,'fi/2m)  is 
the  regular  perturbation  expansion 


r 1/2 

Hence,  the  parameter  i = (b,^)  measures  the  rate 

of  reaction  probability  as  a function  of  the  physical 


T = g 


Zae 

Er^ 


k(r,0,h/2m)  = 5^k"(r,9)  (fi/2m)". 


(4.18) 


n 


and  F(Z)  is  a special  function  that  satifies  (8) 


d^F  _ 7 dF  . 


(4.19) 
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Proceeding  as  above,  we  find  that  and  satisfy  (4.8),  (4.9) 


and  (4.10)  respectively  and  that  k satisfies 
ur,  0 ^ . e,  0 g°  r,  2 , 21  . 


(4.20) 


.0 


The  constructions  given  for  ij/(r,d)  and  f above  can  be  used 
here.  The  value  of  g^,  a constant,  is  determined  by  evaluating 
(4.20)  on  the  limit  cycle: 


0 

2_ 

2 


(4.21) 


i . e . 

g°  = 2/ip^(r^).  (4.22) 

Once  g^  is  known,  k^  can  be  calculated  by  the  method  of  charac- 
teristics (9),  with  initial  data  k^(r,6)  =0  on  r = r^.  If 
4)(r,9)  = on  r = r^,  we  set 

yjh/2m)  = F'(Wyfh7^)=  0.  (4.23) 

when  integrating  (4.19).  With  this  choice  of  k^  and  tp , 

T(r,0)  =0  on  r = r^.  Higher  order  terms  can  be  treated  in  an 
analogous  fashion. 
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AD  ^80  592 
PP  123 

Petersen,  Charles  C '•'The  Sciviet  PortT  tearing  f)pera 
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Ushed  in  the  condensed  version  of  The  Journal  of  the 
AcoustK  Society,  Vol  60,  No  6 Dec  19‘’6  pp  1314 
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and  CylxfTiciits,  Vol  5.  No  1,  May  1975) 

PP  119 

Morgan.  Wilium  P . Jr  . "Beyond  l-olklore  and  babies 
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AD  A 020  443 

PP  143 

Horowitz.  Stanley  and  Sherman.  Allan  (L(  dr..  DSN), 
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AD  A028  189 

PP  148 

McConnell.  James  M . "Military-PutitKal  Tasks  of  the 
Soviet  Navy  in  War  and  Peace."  62  pp. , Dec  1975 
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pp..  Aug  1976  (Reprmicd  frnni  f S \jtj/  InsntuJr 
fV/KreJinys.  Jun  1976)  AD  AO.104‘6 

PP16: 

Fiiedlieim.  Robert  L , "Parliamentary  Diplomacy  " 
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,39th  MORS  Working  (.roup  on  Manjwiwei  and 
Personnel  Planning  Apma)>«'li'  Md  2H  30  June 
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